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We show, via numerical simulations, that the fidelity decay behavior of quasi-integrable systems 
is strongly dependent on the location of the initial coherent state with respect to the underlying 
classical phase space. In parallel to classical fidelity, the quantum fidelity generally exhibits Gaussian 
decay when the perturbation affects the frequency of periodic phase space orbits and power-law 
decay when the perturbation changes the shape of the orbits. For both behaviors the decay rate 
also depends on initial state location. The spectrum of the initial states in the eigenbasis of the 
system refiects the different fidelity decay behaviors. In addition, states with initial Gaussian decay 
exhibit a stage of exponential decay for strong perturbations. This elicits a surprising phenomenon: 
a strong perturbation can induce a higher fidelity than a weak perturbation of the same type. 



PACS numbers: 05.45.Pq, 03.65.Yz 

Manifestations of chaos and complexity in the quan- 
tum realm have been widely explored in connection with 
the correspondence principle between classical and quan- 
tum mechanics ,1] . An example is a system's response to 
small perturbations of its Hamiltonian. Peres '3] con- 
jectured that this response serves as an indicator of chaos 
applicable to both the classical and quantum realms. 
That is, in both realms the behavior of fidelity between a 
state evolved under perturbed and unperturbed dynam- 
ics depends on whether or not the dynamics is chaotic. 

Peres' conjecture found an experimental venue in nu- 
clear magnetic resonance (NMR) polarization echoes. 
In these experiments the initial state of the system is 
evolved forward under its internal dipolar Hamiltonian 
and then inverted by a sequence of radio- frequency pulses 
The inverted Hamiltonian, however, will be not 
be an exact reversal of the internal Hamiltonian due to 
pulse imperfections and interactions with the environ- 
ment. These perturbations reduce the subsequent echo 
amplitude which is the measure of fidelity. 

The polarization echo as a means of studying dynam- 
ical irreversibility was applied in Ref. |^, where it was 
noted that the echo decay behavior as a function of 
time can be exponential or Gaussian, depending on the 
molecule under investigation. The connection between 
these results and the exponential fidelity decay predicted 
for systems exhibiting quantum chaos Q was made in 
Ref. (3. 

Encouraged by these experimental investigations, Jal- 
abert and Pastawski applied semi-classical analysis to 
the evolution of what they termed the Loschmidt echo, 
or fidelity decay. Their analysis showed that for chaotic 
systems, when the perturbation is strong enough such 
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that perturbation theory fails, the fidelity decay is com- 
prised of two exponentially decaying terms. The first of 
these terms is dominant for small errors and can be de- 
scribed by the Fermi golden rule ^ . The second term 
is dominant for strong errors, independent of perturba- 
tion strength, and decays at a rate given by the analogous 
classical system's Lyapunov exponent. 

The identification of a classically chaotic signature in 
quantum systems has led to detailed studies of fidelity 
decay behavior. For quantum systems that are analogs 
of classically chaotic systems, a number of regimes have 
been identified based on perturbation strength. For weak 
perturbations, such that perturbation theory is valid, 
the fidelity decay is Gaussian 0, |^ For stronger 

perturbations, in the Fermi golden rule regime, the fi- 
delity decays exponentially with a rate determined by 
the perturbation Hamiltonian and perturbation strength 
Q> H> 112j IM- I'^ many systems the rate of the 
exponential increases as the square of the perturbation 
strength '3| (see 1^ for an exceptional case) until satu- 
rating at the underlying classical systems' Lyapunov ex- 
ponent 0, or at the bandwidth of the Hamiltonian 01 • 
The crossover between the various regimes [isl EH ll^ 
and the fidelity saturation level 0, ^3 have also been 
explored. Quantum fidelity decay simulations have also 
been carried out in weakly chaotic systems [i^l, and at 
the edge of quantum chaos |2l| . 

Relationships between fidelity decay behavior and 
other quantum phenomena are also found in the lit- 
erature. These include the Fourier transform relation 
between fidelity decay and the local density of states 
Hi nil I issues of classical-quantum correspondence [i^ , 
reversibility [31, and decoherence 25]. We also note that 
fidelity decay studies can be carried out on a quantum 
computer \\A l2fil | , and that the fidelity has been exper- 
imentally determined for a three-qubit quantum baker's 
map on a NMR quantum information processor |27|. 

Studies of fidelity decay in quantum systems have 
spurred interest in the fidelity decay of classical systems 



2 



[28j . For chaotic classical systems it has been shown that 
the asymptotic fidelity decay can be either exponential or 
power-law, analogous to the asymptotic decay of corre- 
lation functions |29j. Faster than Lyapunov exponential 
decays have also been identified [s^] . 

The fidelity decay behavior of quantum analogs of non- 
chaotic or quasi-integrable classical systems has received 
less attention [l3,|3i|33l then its chaotic coun terp art and 
has been the subject of some controversy HE 13- Using 
semi-classical arguments, Prosen |3| demonstrated 
the counter-intuitive result that quantum fidelity decay 
of regular, non-chaotic, evolution is Gaussian, faster than 
the exponential decay of chaotic systems. This was chal- 
lenged by further semi-classical arguments [s^ which in- 
dicated a power-law decay. A proposed resolution j^^l 
differentiates between individual minimum uncertainty 
states, which generally exhibit a Gaussian decay, and av- 
erages over many such states, which may be biased by 
specific states exhibiting power-law fidelity decay behav- 
ior. 

In this work, we explore what causes a quantum state 
undergoing regular quantum evolution to exhibit Gaus- 
sian or power-law fidelity decay behavior. We present 
numerical results demonstrating that the behavior de- 
pends on the reaction of the underlying classical phase 
space to the applied perturbation. Building off classical 
fidelity decay results [s^j we chart the regions of phase 
space containing states with initial Gaussian or power- 
law decay. Within the two regions we show that the exact 
rate of the Gaussian or power-law decay is also a func- 
tion of the coherent state position. In addition, a con- 
nection is presented between fidelity decay behavior and 
the spectrum of the initial state in the eigenbasis of the 
system. Finally, we probe the dependence of the initial 
decay behavior on perturbation strength, Hilbert space 
dimension, and note that, for strong perturbations, there 
exists a transitional exponential fidelity decay behavior 
after initial Gaussian decay and before fidelity saturation. 

Perturbing classical Hamiltonian evolution can affect 
phase space orbits in two general ways: the perturbation 
may distort the shape of the orbit or change the frequency 
of the orbit. Benenti, Casati, and Veble (BCV) pro- 
posed that in the limit of weak perturbations the classical 
fidelity decay behavior is solely determined by the dom- 
inant perturbation effect on the phase space orbits. If 
the dominant effect on a specific orbit is to change its 
frequency, initial wave packets centered in the region ex- 
hibit Gaussian decay (assuming Gaussian wave packets) . 
This is what would be expected from the fidelity of two 
Gaussian wave packets moving in antiparallel directions, 
or at different speeds, along a specific path. If, however, 
the effect of the perturbation is to change the shape of 
the KAM torus, states centered in the region will exhibit 
power-law fidelity decay. BCV note that they expect 
similar results in the quantum realm. 

Here, we provide numerical evidence that the corre- 



spondence between the perturbation's effect on phase 
space and fidelity decay behavior extends to quantum 
systems. Specifically, we show that quantum fidelity de- 
cay behavior depends on whether an initial coherent state 
is centered on a phase space orbit whose frequency is 
changed due to the perturbation, in which case the de- 
cay will be Gaussian, or an orbit whose shape is dis- 
torted by the perturbation, in which case the decay will 
be power-law. Fidelity decay simulations under quantum 
kicked rotor evolution support a suspicion of Ref. [35j |. 
that quantum states are more prone to Gaussian decay 
due to the quantization of the phase space tori. 

The quantum fidelity decay of an initial state \Tpi) is 
given by 

F{t)^mu-%\A)\' (1) 

where U is the unperturbed evolution, Up = C/e^*'^^ is 
the perturbed evolution, S is the perturbation strength, 
and V is the perturbation Hamiltonian. Our numerical 
work is centered around kicked maps with kick strength 
k determining whether the evolution is chaotic or reg- 
ular. For the perturbed evolution we employ the same 
map with a slightly different kick strength. Thus, the 
unperturbed operator is U = U{k), and the perturbed 
operator is Up = U{k + 5k), with perturbation strength 
Sk- 

We begin our study of fidelity decay with the quantum 
kicked top (QKT) a system used in many previous 
studies of quantum chaos iri general and fidelity decay 
in particular H S [H El IHI Ei ■ The classical kicked 
top describes dynamics on the surface of a sphere 

Xt+l = Zt 

Ut+i = xt sin(fcTz) -I- ut cos{kTz) 

zt+i = ~xtCOs{kTz) + ytsin(kTz), (2) 

where fc^ is the kick strength. We choose a kick strength 
oi kx = 1.1 corresponding to quasi-integrable dynam- 
ics. The (j) — 6 phase space of the classical kicked top 
is shown in Fig. ^ where 9t = arccos(2:t) and (j)t = 
arctan(2/t/a;(). The phase space has a stable fixed point 
at {(p = —7t/2,9 = tt/2) surrounded by KAM tori, and 
rotational KAM tori at the S-edges. Another stable fixed 
point is found at {(p = 0,9 ~ Tr/2) encircled by a smaller 
region of stable KAM tori. 

Fig. n illustrates the effect of changing the kick 
strength on the classical kicked top phase space by plot- 
ting orbits of two different perturbation strengths. The 
shapes of the rotational orbits in the regions at the 9- 
edges of phase space and of the tori around the central 
fixed point change significantly while those around the 
fixed point at {(j> = —n/2,9 — tt/2) do not. If corre- 
spondence holds between classical and quantum fidelity 
decay, this observation should alert us as to the likely 
fidelity behavior of coherent quantum states centered in 
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these phase space regions. A power-law decay is expected 
for states centered in the former regions, and a Gaussian 
decay for those centered in the latter region. 

The quantum kicked top (QKT) is defined by the 
Floquet operator 

Uqkt = e-'^^y/^e-^'^^^"'^^. (3) 

where J is the angular momentum of the top and J 
are the irreducible angular momentum operators. The 
Hilbert space dimension of the top is TV = 2J + 1. The 
representation is such that Jz is diagonal. As initial 
states we use minimum uncertainty angular momentum 
coherent states centered around ^i) l3 ^'^^ employ a 
QKT of J = 500 unless otherwise noted. 

For convenience we number the states assuming a 10 
by 10 grid evenly spaced in the 9 and directions as seen 
in Fig. n The lines of the grid are numbered such that 
the number of a state, centered at an intersection of the 
grid, is determined by adding the numerical values of the 
horizontal (numbers on left) and vertical (numbers on 
the bottom) lines. State 1 is thus located at (0 — — 7r/2, 
e = tt/IO) and state 100 at {(j) = 27r/5, 9 ^ tt). In this 
way, the fixed point at {(f) = —tt/2,9 = 7r/2) is number 
41 while the fixed point at {</) = 0,9 = tt/2) is number 
46. 

Fig. El shows that state 52, centered in the region of 
phase space surrounding the fixed point {(f) = — 7r/2, 9 = 
7r/2), exhibits the expected Gaussian behavior (trian- 
gles), and that state 66, near the fixed point (cj) = 0,9 — 
7r/2) (solid line), exhibits the expected power-law decay. 
These states parallel the expected classical fidelity decay 
behavior. 

The existence of states with power-law decay supports 
the contentions of Ref. [sJl in explaining contradictory 
results in regular system fidelity decay behavior. These 
states bias the average to look like a power law which 
is slower than the exponential fidelity decay of chaotic 
systems. Many states, however, exhibit Gaussian fidelity 
decay which may, under certain circumstances, be faster 
than the decay of the corresponding chaotic fidelity. 

The parallel between the fidelity of quantum states and 
their classical counterparts is not, however, the case in 
general. Some states initially centered in areas of ap- 
parent phase space orbit distortion exhibit the Gaussian 
fidelity decay expected from perturbations of an orbit's 
frequency. We show this in the quantum version of the 
system explored classically by BCV IsSl, the kicked ro- 
tor. The classical dynamics of the kicked rotor is given 
by 

Pt+i = pt + (A:fl/27r)sin(27rgt) 

qt+i = qt+Pt+i (4) 

where kn is the rotor kick strength and —1/2 < q,p < 
1/2. For kfi = .3 the phase space of the kicked rotor has 



a stable fixed point at (g = —.5,p = 0) (x) and an unsta- 
ble fixed point at {q = 0,p = 0) (-I-). The phase space is 
divided into two distinct regions, orbits around the sta- 
ble fixed point, and rotational motion, as shown in Fig.l^l 
The orbit at the border between these regions is the sep- 
aratrix. As with the QKT we plot phase space orbits of 
different fc^ to demonstrate the effect of a change of kick 
strength perturbation on different parts of the classical 
phase space. The shapes of certain KAM tori, such as 
the ones just outside the separatrix, exhibit large defor- 
mations while others, such as the inner circles within the 
separatrix, do not. As shown in 35], states in the for- 
mer region exhibit Gaussian classical fidelity decay while 
those in the latter region exhibit power-law decay. 

To study quantum fidelity decay of the kicked rotor, we 
use the unitary operator describing the quantum kicked 
rotor (QKR) 

Uqkr = e-*p''^^e-^'=«™"(2''«)^/'', (5) 

where N is the Hilbert space dimension. For our sim- 
ulations we use QKRs of = 500, 1000, fc/j = .3, cor- 
responding to a classical kicked rotor with non-chaotic 
dynamics, and perturbation strengths Sr — .002, .0014. 
As initial states we use the minimum uncertainty coher- 
ent states described in centered around {qi,pi). 

Fig. 2] shows the fidelity decay of quantum coherent 
states centered in different regions of phase space. One 
(squares) in a phase space region where the effect of 
5ji is primarily to change the frequency of the KAM 
ior\,{q = —A,p = .1) (square in Fig. O, and another 
(solid line) in a phase space region where S]:^ primarily 
distorts the shape of the KAM tori {q = — .l,p = .1), 
(diamond in Fig. ISJ. The former exhibits a Gaussian fi- 
delity decay while the latter exhibits a decay which is 
non-Gaussian and resembles (but is not quite) a power- 
law. An exact power-law does not emerge for any of 
the states simulated for the QKR even for perturbation 
strengths as low as Sr = .0005. Fig.0]also shows a state, 
(q — —.l,p — — A), centered on a rotational orbit outside 
the separatrix which, despite the classical prediction of a 
power-law fidelity decay, exhibits Gaussian decay (trian- 
gle in Fig. Ell. 

As mentioned above, BCV HI] predict a lack of cor- 
respondence between classical and quantum fidelity de- 
cay noting that the quantization of KAM tori tends to 
suppress transitions between them. This would enforce a 
change of frequency as the primary effect of the perturba- 
tion. The lack of an actual power-law decay at any point 
in the QKR orbits implies such a suppression throughout 
the region of rotational motion. Lack of correspondence 
to classical dynamics for rotational QKR orbits has also 
been noted with respect to fidelity recurrences [131 ■ Fi- 
delity recurrences, as seen in Fig. ^ are predicted classi- 
cally in all regions of phase space. Yet, in the quantum 
realm, they do not occur in the QKR rotational tori. 
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FIG. 1: (Color online) Twenty-five classical orbits on the phase space of the classical kicked top, /ct ~ 1-1 (dark), and /ct ~ 1.3 
(light). The same initial points are used to map the orbits so as to demonstrate the efi'ect of the St perturbation. The change 
of kick strength primarily affects the frequency of the orbits around the stable point at {(p = —tt/2, = 7r/2) (x). Thus, under 
a St perturbation, we expect coherent states placed there to exhibit a Gaussian fidelity decay The change of kick strength does 
affect the shape of the central tori around the {(p — 0,9 = n/2) (x) fixed point. Thus, under the St perturbation we expect 
coherent states centered in that region to exhibit a power-law fidelity decay. The same holds for the rotational tori at the edges 
of the phase space. The different markers delineate the fidelity decay behavior of coherent states centered at those points. The 
states marked by circles, up, down, left, and right pointing triangles, five-pointed stars, and six-pointed stars mark states that 
exhibit Gaussian decay behavior at different rates. While we offer no clear, a priori, determination of the decay rate, states 
centered on the same orbit tend to decay at similar rates. Though different orbits may also give rise to similar decay rates, the 
general trend is towards a slower Gaussian as the states move further from the {(p = — 7r/2,S — it/2) fixed point. Power-law 
fidehty decay rates depend on distance from areas with large tori distortions, the region surrounding the central fixed point 
and the dip in the rotational torus at the top and bottom of the figure. The closer a state is to these regions, the slower the 
power-law rate. Gaussian and power-law fidelity decays are shown in Fig. |5|where the shapes correspond to the markers used 
in this figure. 



Ref. [321 notes that these regions have a high density of 
states and the enhanced quantum interference may cause 
the lack of classical correspondence. A similar argument 
can be proposed here, as enhanced quantum interference 
may be especially sensitive to perturbations and cause 
the fast Gaussian fidelity decay. 

In attempt to further understand what differentiates 
states that exhibit Gaussian fidelity decay from those 
that exhibit power-law decay we look at the spectra of 
the initial coherent states with respect to the QKT eigen- 
basis as a function of the extent of the eigenstates in . 
Following Peres Q the extent of a state with respect to 
the operator is defined as 

A|j,i = ^/wIWTWjmW- (6) 

The extent is basically the first term in the power se- 
ries expansion of the fidelity and thus we expect it 
to provide insight into the expected fidelity decay behav- 
ior. We calculate the extent of all the eigenstates of the 



QKT and see how much each of these states contributes 
to a given coherent state. The contribution is quantified 
by an amplitude Aj = \{ipi\<t'j)\'^ where ipt is the initial 
coherent state and (pj is the jth QKT eigenstate. 

On the extremes, the coherent state centered at the 
stable fixed point {(p — — 7r/2, 9 = 7r/2) is primarily {A — 
.95) composed of one eigenstate with A| J^j — 17.3. The 
primary contributors to the coherent state centered at the 
stable fixed point at {(f) = 0,9 — tt/2) are 4 eigenstates 
each with amplitudes of .21 and A| J^j = 353.7. In both 
cases the fidelity barely decays as the state lives in a 
constricted Hilbert space [j, Um ■ 

Most coherent states, however, have significant con- 
tributions from many different eigenstates. In Fig.|Slthe 
contribution to coherent states 52-56 is plotted versus the 
extent of the eigenstates and in the inset the same is plot- 
ted for states 64-67. The general pattern emerging from 
the figure (and from states not shown) is clear. Coherent 
states exhibiting Gaussian fidelity decay have a Gaus- 




FIG. 2: (Color online) Fidelity decay for quantum kicked 
top (QKT) with kr = 1.1, 5t = .001, and J = 500, of two 
different coherent states. One (triangles) placed in the re- 
gion of stable KAM tori surrounding the fixed point (0 = 
— 7r/2, 9 = 7r/2) (state 52, up triangle in Fig.0 every 20 steps 
shown) where the perturbation effects a change in phase space 
orbit frequency, and the other (solid line) placed in the re- 
gion surrounding the {(j) = 0,0 = n/2) stable point (state 
56, diamond in Fig. where the perturbation causes a dis- 
tortion of the phase space orbit. As expected, the former 
displays a Gaussian fidelity decay (dashed line), e~^'^^ , with 
Fg — 4.5 X 10~^, while the latter exhibits a power law de- 
cay cpt~"^ (dash-dotted line) with ap = 1.15 and cp = 950. 
The insets show the change of fidelity decay rate as a function 
of perturbation strength, 6t- The top inset plots St versus 
Cp for states that decay as oc {+), t~^'^^ (diamond), and 
^-1.3 j-*^ j^g^g explained later, different coherent states decay 
with different power-laws). The points on the loglog plot are 
well fit by Cp = .855^^^ (dotted line), .355^^'^^ (dashed line) 
and .855^^'^ (dash-dotted line). From these and numerical 
results of other states we assume the following form for power 
law decay Fp{t) — cp{5Tt)~°'^ ■ The lower inset shows the 
change of the Gaussian rate, Fg, as a function of St (marked 
by shapes corresponding to those in Fig.0. As addressed be- 
low, coherent states centered in different areas of phase space 
exhibit different Gaussian decay rates, with the rate slowing 
as the states move further from the stable point. The depen- 
dence of the rate on perturbation strength, however, is always 
Fg (X St (dotted line). 

sian spectrum of contributions from eigenstates with low 
to middle range extents in J^. As the coherent states 
move away from the low-extent (0 = = 7r/2) 

fixed point (with distance determined by the number of 
passing trajectories) more and more states contribute and 
with lower amplitudes, higher extent, and a more local- 
ized extent range. The Gaussian shape remains until the 
coherent states enter the region of power-law fidelity de- 
cay. States with power-law fidelity decay have a very 
narrow spectra at high extent. The difference between 
the types of states is clearly seen in Fig. |31 and the tran- 
sition between the types of states can be seen in Fig. 
We suggest that the relationship between fidelity decay 
and the extent spectrum may be understood as follows. 




FIG. 3: (Color online) Thirty orbits on the phase space of the 
classical kicked rotors, kp — .3 (dark), and kp = .35 (light). 
The same initial points are used to plot the orbits of both 
maps in order to highlight the effect of a 5p perturbation. 
Well within the separatrix, the effect of the 5p perturbation 
is generally to change the frequency of the KAM tori. Out- 
side the separatrix, however, the 5p perturbation essentially 
changes the shape of the tori. Also shown are the stable, 
(g = — .5,p = 0) (x), and unstable, (g = 0,p = 0) (+), fixed 
points and states used in the text to demonstrate Gaussian 
(square) and non-Gaussian (circle, diamond) fidelity decay. 
Classically, states outside the separatrix are expected to ex- 
hibit power-law fidelity decay. In the quantum realm, how- 
ever, we find that many such states exhibit Gaussian decay, 
an example of which is marked by the triangle. 

States exhibiting power-law decay have a large extent 
in the direction of the perturbation, Jz- When the co- 
herent state is perturbed these states cannot spread out 
much more in the perturbation direction. Rather, they 
interfere with each other and the decay is slow. Coherent 
states exhibiting Gaussian decay, however, are spread out 
in extent space. The perturbation affects each of these 
states differently spreading them out in Jz and causing 
a ballistic decay. As the coherent states move away from 
the low extent fixed point the average extent of these 
states grows and the Gaussian fidelity decay gets slower 
until the transition to power-law decay. This description 
holds for states in the Gaussian and power-law fidelity 
decay regions for the QKT phase space. States at the 
border between these regions and states very close to the 
fixed points have different extent spectra and, thus, ex- 
hibit fidelity decay behavior that is neither Gaussian nor 
power-law. These regions will be discussed below. 

A full exploration of the extent and its relation to fi- 
delity decay is beyond the scope of this paper. However, 
looking at the spectrum of a coherent state as a func- 
tion of the extent of the contributing basis states with 
respect to the perturbation operator, V (or some func- 
tion thereof), can help identify the regions of different 
decay behaviors. 

We now embark on a more extensive exploration of co- 
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FIG. 4: (Color online) Log-log and log-linear plots of three 
coherent state fidelity decay behaviors under evolution of 
the quantum kicked rotor with ku = .3, Sr — .002, and 
A'' — 500. One of the states is centered inside the sepa- 
ratrix at {q = —A,p = .1) (squares every 50 time steps), 
where there is no apparent distortion of KAM tori. The loca- 
tion of this state on the classical phase space is marked by a 
square in Fig. |31 Another state is centered outside the sepa- 
ratrix (q = —.l,p — .1) (solid line, diamond in Fig. |^ where 
there is noticeable KAM tori distortion due to the perturba- 
tion. The former displays a Gaussian fidelity decay (dashed 
line), e~ ooo^* and exhibits fidelity recurrences every 2300 
time steps Is^ . The latter exhibits a decay which is non- 
Gaussian. A power-law oc (dash-dotted line) is plotted 
for comparison. The third state (triangle) is centered outside 
the separatrix (triangle in Fig. |3J but nonetheless exhibits a 

Gaussian decay F{t) = e"^*^^'^ * . This state does not ex- 
hibit fidelity decay recurrences. 



herent state fidelity decay behavior in the regular regime 
of the QKT. To this aim, we have calculated the fidelity 
decay for coherent states spaced throughout the classical 
phase space for a number of perturbation strengths and 
Hilbert space dimensions. A large variety of behaviors 
exist, though we concentrate only on the initial decay 
before any fidelity recurrences. In attempt to organize 
the data in a straightforward fashion we explore fidelity 
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FIG. 5: (Color online) Contribution of QKT eigenstates for 
coherent states 52-56 (left to right) versus extent of those 
eigenstates in Jz with shapes as in Fig. As the states 
move further from the low extent {(fi = —it/2, 6 = n/2) fixed 
point more eigenstates have significant contributions and the 
height of the Gaussian contribution curve decreases. Conse- 
quently, the Gaussian fidelity decay rate decreases. State 56, 
made of only high-extent states, does not have a Gaussian 
extent spectrum and exhibits power-law fidelity decay. The 
inset shows coherent states 64-67 (left to right) which follow 
a similar pattern. 



decay as it relates to the following variables: perturba- 
tion strength, Hilbert space dimension, and position of 
the initial state with respect to the underlying classical 
phase space. We also study hitherto unobserved expo- 
nential decay which may occur after an initial Gaussian 
decay and explore how this decay regime behaves with 
respect to the above variables. We note that an exten- 
sive semi-classical treatment of regular fidelity decay has 
been done in Ref. . Our purpose here is to outline an 
approach based on knowledge of the system's classical 
phase space. 

We first address the dependence of the fidelity decay 
rate as a function of perturbation strength. For coher- 
ent states exhibiting a Gaussian fidelity decay. Fait) — 
g-Tc* ^ numerical simulations verify Tg (x S'^, as derived 
in This dependence is demonstrated in the lower 

inset of Fig. El for 5t = .0001, .0005, .001, .005, .01 and 
J — 500. For coherent states exhibiting a power law 
decay, Fp{t) = cpt~°''', numerical results for the above 
perturbation strengths suggest that cp is cx 5^"^ , from 
which we conclude Fp{t) — c{6Tt)^°'^ ■ The upper inset 
in Fig. 12 demonstrates this behavior with states whose 
power-law decay rate is ap = 1, 1.15, and 1.3. 

To address the fidelity decay behavior as a function 
of Hilbert space dimension we choose one perturbation 
strength, 6t = .005, for J = 100,200,300,400,500. The 
fidelity decay is calculated for coherent states of appropri- 
ate dimension centered at specified points in phase space. 
For states centered in regions of Gaussian fidelity decay 
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FIG. 6: (Color online) Coherent state 56 (diamond in Fig.Ql 
power-law fidelity decay (every 16tli time step plotted) for the 
QKT with kr = 1.1, St = .005, and J = 500 (circles), 300 
(squares), and 100 (dots). Lowering J causes the fidelity de- 
cay behavior to transition from power law to Gaussian. The 
inset shows the rate of Gaussian fidelity decay, To, as a func- 
tion of J for some coherent states. We find a linear relation- 
ship Fg = cgJ with CG = 2.3 x 10"'^, 5.5 x 10"'^, 6.7 x 10"^ 
and 1 x 10~® (bottom to top) shown. 

we find a linear relation between J and Tq, as shown in 
the inset of Fig. El the slope of which depends on the 
coherent state's location in phase space. 

We can thus write the following equation for the Gaus- 
sian fidelity decay behavior: 

FG{t)=e-^^'\ Tg^IgJS^ (7) 

where 7g depends only on the initial coherent state's 
location in phase space and is the only term not a priori 
calculable from our analysis. This result is in consonance 
with the semi-classical approach outlined in [nj . 

For coherent states in regions of power-law fidelity de- 
cay, we find no change in decay rate with J as long as 
the fidelity decay remains a power-law. Thus we write 
the following equation for the power-law decay behavior 

Fp{t) = ciSri)-''^ (8) 

where c and ap depend on the coherent state's location. 
However, as J is decreased the fidelity decay behavior 
does change; it shifts from power-law to Gaussian, as 
shown in Fig. El (for state 56, marked diamond in Fig.^. 
This shift is due to the increasing size of the initial co- 
herent state making it more likely that the states will 
overlap with KAM tori on whom the perturbation effects 
a change of frequency. Coherent states centered more 
deeply in the power law decay region (such as state 77 
marked plus in Fig. ^ have a slower transition to Gaus- 
sian decay when decreasing J. 

As we have seen, the fidelity decay behavior in general, 
and the rate of Fq (t) and Fp (t) specifically, are depen- 
dent on the exact location of the initial coherent state 



with respect to the underlying classical phase space. This 
dependence is emphasized in Fig. ^ by using different 
shapes to mark the center of coherent states exhibiting 
Gaussian (circle, square, up, down, left, and right trian- 
gles, five-pointed star, and six-pointed star) and power- 
law (diamond, dot, *) fidelity decay. States with prac- 
tically equivalent decay rates are represented by the same 
shape, with the rates themselves shown in Fig. [3 Based 
on our simulations we cannot formulate clear cut rules 
for the decay rate of a given coherent state. However, 
we make two observations. First, states along the same 
phase space orbit tend to have similar decay rates. Sec- 
ond, Fg, the rate of Gaussian decay, decreases as the 
states get further from the (</> = —tt/2,9 = tt/2) fixed 
point (with distance measured by the number of trajec- 
tories between the fixed point and the center of the co- 
herent state) . We have already seen consequences of this 
latter observation in the extent spectra. Similarly, for 
states exhibiting power-law decay behavior, the power 
increases as the states get further from regions of large 
KAM torus distortion. 

We now explore the fidelity decay of states at the bor- 
der between the Gaussian and power-law phase-space re- 
gions and of states close to the fixed points. These regions 
exhibit a variety of decay behaviors which are reflected 
in the extent spectrum. Looking at the transition from 
Gaussian to power-law decay away from the fixed points 
we note that the transition is a rather smooth one in 
both the decay behavior itself and the extent of the con- 
tributing eigenstates. Fig. [3 displays these behaviors for 
coherent states 9 = Att/5 and (j) ranging from to tt/IO 
(between states 76 and 77 of Fig. For the fidelity, 
the decay slows as the state leaves the region where the 
dominant perturbation effect is on the frequency of the 
orbits and enters the region where the dominant pertur- 
bation effect is on the shape of the orbits. In the extent 
spectra the transition is manifest by the Gaussian shape 
narrowing on the side of large extent eventually becoming 
almost flat except for a tail reaching towards the higher 
extent eigenstates. Similar behavior is found for other 
states in the border region between Gaussian and power- 
law fidelity decay. 

The region surrounding the {4> = 0,0 = 7r/2) fixed 
point contains states exhibiting fidelity decay behaviors 
not described by a Gaussian or power-law and extent 
spectra different from those seen above. Starting with 
the coherent state centered at the fixed point, the fidelity 
oscillates close to one as the coherent state is comprised 
almost entirely of the highest extent QKT eigenstates. 
As the coherent states move away from the fixed point the 
highest extent eigenstates still gives the largest contribu- 
tions while the next highest extent eigenstates give in- 
creased contributions. The fidelity decay in these regions 
starts off as a power-law but exhibits a second stage of 
Gaussian decay similar to edge of quantum chaos decays 
[2l| . Moving further, eigenstates with lower and lower 
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FIG. 7; (Color online) Fidelity decay of coherent states with 
6 = 47r/5 and (f) ranging from (state 76) and tt/IO (state 
77) for the QKT with kr = 1.1, St = .001. The states 
shown are for <f) = vr/lOO, Svr/lOO, Stt/IOO, 67r/100, 77r/100, 
and 97r/100 (bottom to top). The decay can be seen to transi- 
tion smoothly from Gaussian to power-law. The inset shows 
the extent spectrum of the states at = Svr/lOO (circles), 
Gtt/IOO (diamonds), and Qvr/lOO (x). As the power- law fi- 
delity decay region is approached, the Gaussian spectrum gets 
filled in and narrows on the side of higher extent becoming 
almost fiat except for a tail of high-extent states. 



extent become dominant. However, the extent spectrum 
is not Gaussian, as would be expected for states exhibit- 
ing Gaussian fidelity decay, but is extended on the side 
of lower extent eigenstates. At this increased distance 
from the fixed point, the first-stage fidelity decay transi- 
tions from power-law to Gaussian, and the second-stage 
decay transitions from Gaussian to power-law decay. Fi- 
nally, the initial Gaussian fidelity decay flattens into one 
stage of power-law decay while the extent spectrum con- 
tinues flattening in the direction of higher extent states 
while forming a complicated flattened bulge at lower ex- 
tent states. All of these behaviors are exhibited in Fig.|51 
Similar behavior is found in regions of coherent states 3-5 
and 87-89. 

Coherent states in the region surrounding the fixed 
point at (0 = —tt/2,0 = 7r/2) exhibit behavior that is 
slightly different from the states in the region surround- 
ing the {(j) — Q.6 — 7r/2) fixed state. At the fixed point 
the fidelity simply oscillates close to 1. As the coherent 
states move away from the fixed point the oscillations be- 
come larger in amplitude, the recurrence time increases, 
the initial decay becomes more Gaussian, and the max- 
imum fidelity reached on the recurrence is lower. This 
continues until full Gaussian decay behavior emerges. 
The extent spectra refiect this behavior, going from a 
dominant low extent state to an eventual Gaussian shape. 

We also note the presence of states with unexpected 
fidelity decay behavior in the QKR. Fig. |^ shows exam- 
ples of = 1000 coherent states under QKR fc/j = .3, 




FIG. 8: (Color online) The top panel shows the fidelity decay 
of coherent states with 6 ranging from 7r/2 to Gvr/lO and = 
for the QKT with kr ~ 1.1, St ~ .001. The states shown are 
for 9 = 7v/2 + 0, 7r/100, 27r/100, Stt/IOO, 47r/100, and Svr/lOO, 
(top to bottom). Close to the fixed point the fidelity oscil- 
lates close to one. As the states move further away from the 
fixed point, the fidelity decays in two stages, starting as a 
power-law and becoming a Gaussian. For comparison we fit 
the initial decay of the 6 = 7r/2 -I- 27r/100 state with a power- 
law (X (dotted line), the — tv/2 + Stt/IOO state initial 
decay with a power-law oc t~^'^ (dashed line), and the sec- 
ond stage decay with a Gaussian e"'"''''* with Fss = 8 x 10~* 
(chained line). Moving further away, the initial power-law 
becomes more Gaussian while the second stage starts fiat- 
tening to power-law. The inset shows every 60 steps of the 
fidelity decay for the states (f} = 0, 9 = 7r/2-|-67r/10 (triangles), 
77r/100 (diamonds), 87r/100 (-I-), 97r/100 (x), and tt/IO (state 
56, squares). Here, the initial decay starts off as Gaussian and 
rebounds into a power-law decay. For comparison we plot the 
Gaussian e~'""^* with Fa = 3.9 x 10~^ (chained line) and the 
power-law oc (solid line) . As the states continue to 

move away from the fixed point, however, the Gaussian flat- 
tens until there is a single-behavior power-law decay, cc t~^'^^ 
(dashed line). The lower panel shows the extent spectrum of 
the states with 6* = 7r/2-|-57r/10 (o), 671/ 10 (triangles), 77r/10 
(diamonds), Ott/IO (x), and state 56 (squares). As the coher- 
ent states are moved further away from the flxed point the 
spectra flatten at higher extent eigenstates and a bulge grows 
at lower extent eigenstates. 
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FIG. 9: (Color online) The upper plot shows two coherent 
states evolved under the — 1000, kR — .3 QKR with 5r — 
.0014 that exhibit initial exponential decay despite the fact 
that the QKR is in the regular regime. The lower plot shows 
three coherent states under the same QKR evolution which 
exhibit initial Gaussian decay that transitions to exponential 
decay. For two of the states the fidelity freezes after the initial 
Gaussian decay. The fidelity of all three states exhibit echo 
resonances similar to those discussed in Ref. |3l| for ergodic 
perturbations. 



5fi = .0014 evolution that exhibit initial exponential fi- 
delity decay (top plot), though the QKR is regular, and 
fidelity 'freeze' as discussed in Ref. [31|, though the per- 
turbation is non-ergodic. 

Beyond the initial Gaussian fidelity decay of some co- 
herent states, there may exist a second, slower, stage 
of exponential fidelity decay behavior, F{t) — CEe~^'^ , 
before saturation. This stage is prevalent for strong per- 
turbations but disappears for smaller perturbations (or 
smaller Hilbert space dimension with the same pertur- 
bation strength) . The specifics of this exponential decay 
depend strongly on the phase space location of the ini- 
tial coherent state. The slower exponential decay of this 
second stage gives rise to an exciting phenomenon: a 
stronger perturbation leading to a higher fidelity than a 
weaker perturbation of the same type. 

The golden rule exponential fidelity decay term men- 
tioned above for chaotic systems exists also in regular 
systems [s^ . We do not identify this term with the expo- 
nential observed here since, as we show, the exponential 
here is strongly dependent on the initial state. 

As with the initial fidelity decay behavior we attempt 
a systematic numerical analysis of the second stage expo- 
nential decay. We first study the exponential as a func- 
tion of perturbation and then explore the effect of the 
location of the initial coherent state. 

Fig, irni demonstrates that weaker perturbations (-f,o) 
exhibit no exponential fidelity decay stage. Rather, the 
Gaussian decay continues until fidelity saturation. As 
the perturbation strengthens the second-stage exponen- 




t 



FIG. 10: (Color online) Coherent state 54 (circle in Fig. 
fidelity decay (every 5th time step plotted) for the QKT with 
kx = 1.1, and J — 500, for perturbation strengths St ~ .01 
(squares), .0075 (diamonds), .005 (triangles), .0025 (+), and 
.001 (o). The fidelity decay of the stronger perturbations 
show a stage of exponential decay after the initial Gaussian. 
This stage is fitted by F{t) = 12e" °^^' (dash-dotted line) for 
St = .01, F{t) = .025e-'"'"5* (solid line) for St = .0075, and 
F{t) = 7e- f'''"'2« (dashed line) for St = .005. For the strongest 
perturbations, St = .01, .0075, there is also a period of trans- 
fer between the two fidelity decay regimes, while for S — .005 
this transition period is non-existent. For weaker perturba- 
tions the exponential stage of fidelity decay disappears alto- 
gether. The faster Gaussian decay of the weaker perturba- 
tions and the wide range of exponential decay rates, leads to 
the counter-intuitive result that a stronger perturbation leads 
to a higher fidelity as seen in the region 200 < t < 600 for 
the St = .01 decay and across the whole plotted region for 
the St ~ .0075 decay. There is no clear correlation between 
perturbation strength and exponential decay rate. 

tial emerges. In addition, there exists a transition pe- 
riod between the two decay behaviors. Thus, weaker 
perturbations lead to longer times of Gaussian fidelity 
decay during which the fidelity of stronger perturbations 
may have already transferred to the slower exponential. 
In this way, there may be a significant amount of time 
in which the fidelity of the stronger perturbation (di- 
amonds, triangles) is actually higher than that of the 
weaker perturbation (squares, +). This phenomenon is 
shown in Fig. [TOI for the kr = 1.1, J = 500, QKT. 
An exponential region of decay is manifest for pertur- 
bation strengths 6t = -01, .0075, and .005, but not for 
weaker perturbations St — .0025 and .001. Thus, for 
times t > 300 the fidelity of at least one of the stronger 
perturbations is higher than the fidelity of a weaker per- 
turbation. 

The rate of the exponential also depends on the pertur- 
bation strength, St- However, our numerical simulations 
do not show any simple relationship between the expo- 
nential rate and the perturbation strength. Rather, the 
rate changes drastically ranging from practically zero, 
decay freeze, to a fast exponential. This also allows 
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a stronger perturbation to have a higher fideHty than 
a weaker one. This is exempUfied in Fig. ^] by the 
S — .0075 perturbation whose fideHty decays very slowly 
and, thus, after a time, is higher than the fidelity of all 
of the weaker perturbations. 

The possibility of a stronger perturbation leading to 
a higher fidelity may have important consequences for 
quantum simulations in which a quantum system is try- 
ing to simulate a given dynamics: a strong error in the 
dynamics may be easier to correct via quantum error 
correction than a weak one. This could allow for an in- 
teresting error-correction scenario. A weak error strongly 
affecting a system should be purposely strengthened so as 
to more accurately perform the desired simulation. This 
would be especially significant in a case where the effect 
of the error is too strong for conventional quantum error 
techniques but can be brought below the error-correction 
threshold if the error is strengthened. 

The existence of the exponential fidelity decay region 
may be related to the quantum freeze of fidelity discussed 
in for ergodic perturbations. In fact, the lower plot 
of Fig. El displays the fidelity decay of coherent states 
evolved by the QKR which actually freeze, though the 
applied perturbation is non-ergodic. 

The region of exponential fidelity decay varies dramat- 
ically with the location of the coherent state on the un- 
derlying classical phase space. This is displayed in Fig. 1111 
where a wide range of exponential decay rates are found 
for different coherent states though they undergo equiv- 
alent evolution. In addition, the time of the transition 
period from Gaussian to exponential varies from state to 
state. 

In conclusion, we have provided a numerical study of 
fidelity decay behavior for coherent states in a quantum 
system whose classical analog is quasi-integrable. We 
find that the initial fidelity decay behavior and rate will 
depend on the perturbation strength, Hilbert space di- 
mension, and initial coherent state location. The quan- 
tum fidelity decay behavior generall y co rresponds to the 
classical fidelity decay explored in and the predic- 
tion therein: quantum states tend more towards Gaus- 
sian decay due to the quantization of the phase space 
orbits. In addition, we show that the spectrum of the 
initial coherent state with respect to the system eigen- 
state extent contains information regarding the fidelity 
decay of that state. Finally, we find that after initial 
Gaussian decay behavior, there may be a second stage of 
exponential decay for strong perturbations. The rate and 
inception of the exponential decay depend on the pertur- 
bation strength and location of the coherent state. The 
existence of this second-stage decay behavior leads to the 
counter-intuitive result that stronger perturbations may 
lead to higher fidelity, a phenomenon which may be im- 
portant for quantum computation. 

The authors acknowledge support from the DARPA 
QuIST (MIPR 02 N699-00) program. Y.S.W. acknowl- 
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FIG. 11: (Color online) Fidelity decay for different coherent 
states (every 5th time step plotted) for the quantum kicked 
top with kr = 1.1, J = 500, and 5t = .01. After the initial 
Gaussian the fidelity decay transitions to an exponential, the 
rate of which depends strongly on the location of the coherent 
state with respect to the underlying classical phase space. 
The states shown are 53 (diamonds), 54 (o), 55 (triangles), 
and 74 (squares). All of these states have initial Gaussian 
fidelity decay, as seen in Fig. |21 which transitions into an 
exponential decay. The current figure exhibits the wide range 
of exponential decay rates and transition times that can occur. 
Starting with the lowest plot we find exponential decays of 
F{t) = e- °55* (diamonds), 10e-°3«t (circles), .175e- °^^' 
(triangles), and .066" °°°'^^' (squares). These exponentials are 
displayed by dashed lines in the figure. 
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